arXiv:1503.08869vl [physics.plasm-ph] 30 Mar 2015 
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The spatial-temporal evolution of the purely transverse current filamentation instability is 
analyzed by deriving a single partial differential equation for the instability and obtaining the 
analytical solutions for the spatially and temporally growing current filament mode. When 
the beam front always encounters fresh plasma, our analysis shows that the instability grows 
spatially from the beam front to the back up to a certain critical beam length; then the 
instability acquires a purely temporal growth. This critical beam length increases linearly 
with time and in the non-relativistic regime it is proportional to the beam velocity. In the 
relativistic regime the critical length is inversely proportional to the cube of the beam Lorentz 
factor 7 o{,. Thus, in the ultra-relativistic regime the instability immediately acquires a purely 
temporal growth all over the beam. The analytical results are in good agreement with 
multidimensional particle-in-cell simulations performed with OSIRIS. Relevance of current 
study to recent and future experiments on fireball beams is also addressed. 

I. INTRODUCTION 

The interaction of energetic particle beams with plasmas is ubiquitous in laboratory and 
in astrophysical scenarios, and so are beam-plasma instabilities such as Weibel [I], current 
filamentation [Hi and two stream The first two instabilities, also referred as Weibel-like 

instabilities, are electromagnetic in nature and arise due to the anisotropy in the momentum 
distribution of the electrons, protons and ions. Specifically, for the current filamentation instability 
(CFI) the role of the velocity anisotropy is played by the counter-streaming flow of the particle 
beams. These instabilities generate exponentially growing magnetic fields, providing one of the 
possible mechanisms for generating near equipartition magnetic fields in extreme astrophysical 
scenarios, such as Gamma Ray Bursts (GRB) |^, and are also closely associated with the 
formation of relativistic Weibel mediated collisionless shocks [7j in space [8] and laboratory 
plasmas [MS]. Recently, the onset of the current filamentation instability was experimentally 
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observed in counterstreaming plasmas in high power laser experiments [nHi3]. Experiments on 
laser wakefield acceleration have also reported the filamentation of the accelerating particle beam 
as it interacts with the background plasma [H]. These instabilities provide an efficient way of 
restoring the isotropy in collisionless plasmas, since the energetic particles scatter off the self 
generated magnetic fields by which the longitudinal momentum is transferred to the transverse 
momentum. 

The available theoretical models for CFI are restricted mainly to a purely temporal analysis 
I21E1 fT5] and do not capture any spatial characteristics of the instabilities, which can be very 
relevant for finite size systems HMSKinHIH]. 

In this paper we obtain the relativistic spatial-temporal solutions for the unstable transverse CFI 
modes in cold plasmas. Our work and approach are inspired by Refs. |191120) . A single differential 
equation is derived to model the instability, considering only the electron response, ignoring the 
finite transverse dimension effects, considering a semi-infinite plasma slab and including the effects 
of a beam density ramp. For a step-like Heaviside beam profile analytical solutions are obtained 
for physically relevant and realistic initial conditions. We further obtain the quasi-static and 
asymptotic behavior of the solutions. The theoretical results are compared with multidimensional 
particle-in-cell (PIC) simulations using OSIRIS [21] . Such spatial-temporal analysis, shown in the 
later part of this paper, is relevant to the jets emitted by the x-ray binaries where the velocities of 
the jets are in relativistic range ~ 0.6c |22| where spatial effects in the CFI modes are significant, 
or to the fireball-like beams [161118] interacting with the plasma. 


II. THEORY 

We consider a two dimensional (2D) slab geometry, where a relativistic beam with velocity 
vohZ and density nQbF{z,t), where F{z,t) is the initial density profile of the beam, is propagating 
in a stationary plasma comprised of cold electrons and immobile ions with homogeneous plasma 
density nop. We analyze the stability of a transverse CFI mode with wavenumber k, and vector 
potential A = A{'ip,T)z exp[ikx], where V' = — z and t = t, which satisfies the Coulomb 

gauge condition V • A = 0 by solving the wave equation (v^ — d^jc^)A = —ATrJfc. Under the 
slow envelope approximation |(9^A| <C \kA\, the governing equation for the vector potential of the 
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electromagnetic wave driven by a current density Jz can be written as 
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where J = —e[nofe-F(V')'yife + nQpVip + nifcUob] is the current density driving the vector potential 
A = Az, 7o6 = \j is the beam Lorentz factor, and c is the velocity of light in vacuum. 
The suffixes 0 and 1 are the 0*^ and 1'^* order perturbed values of the plasma (p) and beam (b) 
parameters defined as plasma electron velocity and density Vp = Tip, Up = nop + nip, and beam 
electron velocity and density Vb = vob + vib, nb = nobF{il!) + nib. The chosen vector potential 
perturbation will generate a magnetic field B = V x A = —ikA{^lJ,T)y exp[ikx] in the y direction. 
Resorting to the fluid equations of motion of a two-species-plasma (plasma electrons and beam 
electrons), using the continuity equation and the equation of momentum conservation for the 
relativistic beam and the stationary background plasma electrons, and restricting to the first order 
values in the weakly coupling limits by ignoring the {dr + VQbd^) term with respect to kvQb, the 
perturbed quantities can be written as. 


Tip =- A, 


me 


drhobVib + lob'^obVibzz] = - ikvobAF^^x - {dr + Vobd^)A, 

me me 


drUib = 


nobeFjij!) 

mc^lob 


^ ^06^ T 2 T Vobd^p') 


Tie 


ikx 


( 2 ) 


Incorporating Eq. Q in Eq. Q by taking the second order r derivative of Eq. Q and further 
neglecting the higher order derivatives dr,d.^dr, d'^ when compared with k'^c^, we obtain, 


[d^, + QFii;)dl,-TlFii,)]A = 0, 


( 3 ) 


where Eq = kvobUpb/V^fobD, D = k'^c^ + + oJ^bhlb^ Q = ‘^^Ib'^ob/{llbD), ujpb = n^b^F/{me^) 
and ujpp = i/nopc'^j (meo). Considering an infinite beam [E('!/>) = 1] and ignoring the second term 
in Eq. Q we retrieve the well known purely temporal evolution of the system with growth rate 
Eq. Interestingly, the equation obtained by Mori et al. |19j (Eq. 10 in Ref. |T9]) to analyze the 
spatial-temporal evolution of Raman forward scattering has the same form as Eq. ([^ obtained here 
for the case of CEL Equation Q can be solved numerically for any general beam profile; however, 
to obtain analytical expression, we assume E(?/)) = LI(V’)) where F[{^) = 0 for ijj < 0, F[{il;) = 1 
for > 0 is the Heaviside function. Respecting causality, we can impose H = 0 for r < 0, and 
define the double Laplace transform of A{t, 'll;) as 
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Thus, by doing the double Laplace transformation of Eq. (§, according to Eq. Q, we obtain the 
field expression in Laplace space as 


A{a,P) = 


QA{0, 0) -ia{l + ^0^ A(0, /3) - iaQA{a, 0) - 1^(0, j3) 


(5) 


+ Qal3 + rg 

where A(0,0), ^(0,/3), A{a,0) and drA{0,f3) are the Laplace transforms of A{t = 0,'i/j = 0), 
A(t = 0, V'), A{t, = 0) and drA{T = 0, respectively, which are the required initial conditions. 
The field A{t, i/j) can be obtained by performing inverse Laplace transformation of A(a, /3), defined 


as 
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where are chosen such that the contour from oo — to “OO — io'(^a,i3) lios below all the 

singularities. For the sake of simplicity we consider the following realistic initial conditions. 


A(t, Ip = 0) = A{t = 0,^p) = An, and drA(T, ip = 0) = 0, 


(7) 


which considers that at r = 0, there is an initial constant noise source throughout the beam and 
for r > 0 the beam front (ip = 0) always encounters fresh plasma, and hence a constant noise 
source. The noise source for most instabilities are considered to be associated with the thermal 
fluctuations, and if thermal fluctuations have no time or space dependence, the constant noise 
source assumption holds correct. Longitudinally modulated or time dependent noise amplitude 
can be some of the forms of noise source that should be considered and the detailed analysis of the 
effect of different noise sources on the CFI spatial-temporal evolution will be addressed elsewhere. 
The above conditions yield ^(0, f3) = An/ (*/?), A(a, 0) = An/ (ia), A(0,0) = An and drA{{), /3) = 0, 
which leads to the solution of Eq. (©, by inverting Eq. ([^, as 
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where Ij is the order modified Bessel function of the first kind [23]. Neglecting the term in 
Eq. Q leads to the quasi-static solutions, which are valid at the beam front for ip <C Qt, as 


A(t,iP) = AriH(r)H(V^)Io 



(9) 


Moreover, and using the stationary phase method which gives the impulse response due to a 
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FIG. 1: Evolution of the field AjAn, (a) comparing (i) the exact [Eq. ([^], (ii) the quasi-static [Eq. Q], 
and (hi) the asymptotic solutions[Eq. @] of Eq. (|^ for the spatial-temporal evolution of purely transverse 
CEI modes at r = 20uj~p , and (b) showing the temporal evolution of Eq. ^ at different times and along 
the beam. The beam is propagating with 70 ;, = 1.25 along the z direction in the equally dense (nob = nop) 
plasma. The wavenumber of the CEI mode is fc = 0.628wpp/c taken to get maximum Q = 0.32c at 70 b = 1.25. 
The arrow pointers on the lines indicates Tsat, which varies as 0.32ct. 

localized initial disturbance at r = 0 and i/j = 0, the asymptotic solution for A(t, ip) at large r can 
be written as m 

Air,Ip) = Anexp(2ro7/rVVQ)- (10) 

The partial differential equation governing the CFI [Eq. ([^] and its exact solution [Eq. ([^] are 
valid for Qt 3 > l/fc, whereas the asymptotic solutions are valid for Qt S> 2ro/(Qfc). 

It is evident from figure [^a), that the quasi-static and asymptotic methods [{ii) and {Hi) in 
figure (a)] fail to capture the spatial saturation of the instability at the back of the beam as 
demonstrated by the full exact solution of Eq. ([^ [line (i) in figure . This specific characteristic 
is also evident in the simulation results to be discussed later in this paper. It is worth mentioning 
here that the asymptotic approach, used extensively for spatial-temporal analysis of the longitudinal 
beam-plasma instabilities [ 21 ] , overestimates the growth and does not seem to give correct spatial 
characteristics for the transverse instabilities, specifically for the CEI discussed in here. 

At the beam front, for ip <C Qt the quasi-static solutions given by Eq. (© match well with 
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FIG. 2: Dependence of cross coupling coefficient Q with beam velocity vob at nob = uop- The dashed line 
has a slope of 1, plotting Q = vob- We note that large Q implies a larger region behind the beam front where 
the spatial-temporal behavior is significant. 


the exact solutions described by Eq. Q. The mildly relativistic (job = 1-25) solutions for A{T,'ip) 
[Eq. Q], presented in figure[^(b) with respect to if: for different times r, indicate that the filaments 
grow spatially from the beam front (ip = 0) to the back until the transition point ipT = Qt. After 
the transition point the instability grows in a purely temporally fashion. We define the beam length 
over which the instability grows spatially as Lgat = Qt [identified as vertical arrows in each line 
in figure (b)]. Beyond this length the instability grows with spatially constant temporal growth 
rate Eq. 

As observed from the previous discussions, and from Eqs. and ( fTo] ) , the spatial-temporal 

behavior depends on the cross coupling coefficient Q. To address this, in figure we analyze the 
dependence of Q on the beam velocity vob for different CEI wavenumbers k at nob = For 
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kc ^ ujpp, LOpb, the maximum value of Q is achieved for vob = 0.5c (70b = 1-15), and varies as 
Qmax - 0.65nob/{nQpk‘^c^). For kc < Upp = ujpb, Qmax - 0.4 at vob = 0 . 6 c ( 70 b = 1.25). In the 
non-relativistic scenario 70 b — 1 and kc <C ojpp = Upb, Q = uob, which is shown as a dashed line in 
figure]^ At higher 70 b 3> 1, A; >> Uppjc or nob 'C riQp, Q tends to 0, and the instability acquires a 
purely temporal behavior. 

We have also considered a beam profile with F{^p) = 1 — exp[—'0^/L^], of direct relevance for the 
comparison with simulations. For such beam configurations the numerical solution of Eq. Q gives 
the same spatial-temporal behavior predicted by Eq. Q but with an enhanced saturation length 
Lgat — Qt + L. The results are compared in figure]^ which will be discussed in connection with the 
simulations performed in the next section. Eor E(^) = 1 — exp[—'0^/L^], the beam density profile, 
and hence the effective temporal growth rate oc ro\/E, reaches the maximum density growth rate 
on the spatial scale length L. Thus, in presence of a density ramp the spatial evolution of CEI can 
be attributed both to the beam density spatial profile, and to the cross coupling term. If L ;§> Qr, 
we may ignore the contribution from the cross coupling term (Qd^^) in Eq. (j^, resulting in a field 
varying as A = Aocosh[FoT\/l — thus determining an extra condition for the relevance 

of the spatial-temporal effects of the current filamentation instability. 

III. SIMULATIONS 

In order to confirm and to explore the theoretical results presented above, we have performed 
2D PIC simulations using OSIRIS |2Ij . We consider a scenario where a relativistic beam is 
propagating through a cold plasma, where the beam and the plasma are both comprised of 
an electron-proton neutral plasma (the temperature of the background plasma is set to zero). 
The simulation box, with dimensions 250 x I00(c/a;pp)^, is divided into 12500 x 5000 cells 
with 3x3 particles per cell per species. The beam has a gaussian density ramp at the front, 
E(V') = I — exp[—where L is the length of the gaussian ramp at the beam front. When 
L —)• 0 this profile mimics a sharp rising beam front with a Heaviside function profile which was 
considered to obtain the analytical solutions [Eq. ([^]. We seed the instability, in order to analyze 
a single CEI mode with wavenumber fegeed) with a small magnetic field perturbation of the form 
By{T = 0) = 5BQCOs{kseedx)-, where 5Bq = 5 x lQ~^mcujpp/e. 

The comparisons between the CEI magnetic field By evolution predicted by the theory [numerical 
solutions of Eq. Q] and the fields observed in the simulations, plotted in figure show that 
the solutions given by Eq. ^ are the most suitable model, among the three models discussed 
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FIG. 3: Comparison between simulations and theory: (a), (b), (c) and (d) represent the magnetic field 
amplitude averaged over transverse dimension of the simulations (solid dark line), theoretical field estimates 
(dashed dark line) and beam profile in the simulations in arbitrary units (red/light solid line) at time 
Tsim = 8.48, 11.30, 14.13, and 16.95uj~p respectively for 70 b = 1.25. To maximize the Q and minimize 
the effect of density ramp on the spatial-temporal evolution of CFl mode we have taken k — 0.628wpp/c, 
nob = n-op and L = Q.'oc/uipp. (e) represents temporal evolution of the Lgat (o) compared with the theoretical 
estimates [Eq. (|^] (dashed line: Q = clrLsat) and (f) represents the logarithmic growth of the field in the 
region of purely temporal growth {ip ^ V't)- Simulation results in o, and theoretical growth rate Fq in dash 
line]. 
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here to predict the spatial-temporal growth of the CFI along the beam, as expected. In the 
simulations the instability starts to grow after a relaxation time (rreiax ~ l-GSUpp for this particular 
simulation) necessary for the self-consistent electromagnetic fields and the electromagnetic noise 
to adjust to the initial flow condition. Thus, for the comparison with the theory, the time is 
re-normalized to r = Tsim — Treiaxj where Tsim is the simulation time. One can observe in figure 
that the theoretical estimate for the CFI magnetic field, given by Eq. ([^, matches well with the 
magnetic field profile observed in the simulations. We analyze the variation of saturation length 
Lsat with time r in figure (e) obtained from the simulations. The rate at which Lgat increases 
with time r is Q = 0.32c, which is equal to the theoretical value of Q = (OrLsat) as predicted by 
our model. Beyond the beam length Lgatj the magnetic field amplitude is spatially constant and 
grows temporally with growth rate Fq, as predicted by the theory. 

The longitudinal modulations, with wavelength ~ Ap, seen in the simulations of figure]^ are due 


100 

80 

^60 

o 

« 40 
20 

0 


BJ- 




—3 X 10 ® By[mcu:pp/e] 3 x 10 ® 




104 106 108 no 108 no 112 114 

z\c/ujpp\ 


FIG. 4: Effect of 70 ;, on spatial evolution of the current filamentation instability: Snap shots of magnetic 
field By at time = ll.SOWj^p^ for (a) yoh = 1.25 and (bjyob = 10, demonstrating that at high Lorentz factor 
the spatial properties of the instabilities are negligible. Other parameters are: fcgeed = 0.628a;pp/c, nob = nop 
and L = O.bc/uipp. The white line shows the beam profile in the simulations. After time r > 
longitudinal modulations on the CFI becomes significantly strong to visualize the spatial saturation of the 
fields. 


to the growth of the longitudinal instability seeded by the sharp rising beam density profile at 
the front. In the simulations, the longitudinal electric field modulation is observed, but confined 
only in the front portion of the beam. The remainder of the beam does not show any sign of the 
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FIG. 5: The current filamentation fields growing along the beam at time Tsim = in the simulations 

performed with varying the size of density ramp L for 7 o(, = 10 and k = G.Oupp/c. The longitudinal mode 
responsible for the modulations at the front of the beam are suppressed for longer density ramps due to the 
reduced noise level for excitation of the longitudinal modes. 


longitudinal electrostatic instability. We attribute this to the fact that, in a similar way as for the 
CFI, the longitudinal instabilities also have a spatial-temporal nature [23]. This also demonstrates 
that a full understanding of this scenario requires the combined analysis of CFI and longitudinal 
electrostatic instabilities. 

The magnetic field snapshots in the x — z plane, shown in figure indicate a sharper rise in the 
magnetic field at the beam front for higher yofe [figure (b)] as compared to lower 706 [figure 
(a)], which further validates the theory, since at high jQb the cross coupling term Q, and thus 
the saturation length Lgat, decreases as ~ l/ 7 ofe for a given time r. Since Fq oc l/y^yob, the field 
amplitude at the back of the beam in figure (b) (high 70 b) is weaker as compared to the field in 
figure]^ (a) ( low job), also as predicted by the theory. 

As observed in figure]^ where the transversely averaged By is compared for various density ramps, 
on increasing the ramp size L, and thus reducing the initial seed for the longitudinal modes, the 
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simulation results show that these longitudinal modulations on the purely transverse CFI modes 
(at the wavelength ~ Ap) can be suppressed. As the time progresses (not shown in the paper) the 
longitudinal modes as well as other faster growing CFI modes start to play an important role and 
their interplay in the nonlinear stage becomes significant. 


IV. DISCUSSION AND CONCLUSIONS 

To summarize, in this paper we have derived a single differential equation modeling the 
spatial-temporal evolution of the purely transverse current filamentation instability. For relevant 
initial conditions exact analytical solutions have been obtained and compared with the analytical 
solutions under the quasi-static and asymptotic approximations. The validity of the model was 
demonstrated by comparing it with 2D PIC simulations in OSIRIS [21]. In a setup of a cold 
relativistic beam propagating in a uniform cold plasma the instability grows from the beam front 
to the back, acquires maximum value at the critical beam length Lsat = Qt at given time r and 
then grows in a purely temporal manner for the rest of the beam length. 

For relativistic fireball electron-positron beams dZl 18 ] undergoing current filamentation in an 
electron-ion plasma, the cross coupling coefficient Q is enhanced by a factor of 2 (Q —)• 2Q) 
and the purely temporal growth rate is enhanced by a factor of V2 (Fq —)■ \/2Fo) due to the 
contribution from the current driven by the velocity and density perturbations of the positrons in 
the beam. However, in the relativistic regime this enhancement is not sufficient to balance the 
1 / 7 q^ dependence of Q on the beam Lorentz factor 70 b. As a result for an ultra-relativistic 29 GeV 
fireball beam m with nofe = nop, Q — 2 x 10 and thus the spatial evolution to the CFI can 
be attributed only to the density gradient scale length. For the recent experiments with a OOMeV 
electron beam |T 6 |, Q ~ 1.2 x 10“®c, thus again suggests only purely temporal growth of the CFI 
is present along the beam. However, in the case of moderately relativistic fireball beams Q can be 
significantly enhanced. For instance, in the case of Sarri et al. US] with 7 ob = 15, nob = lOnop 
and considering a density ramp of L = 0.22cfu}pp, the cross coupling coefficient can be Q ~ 0 . 01 c, 
which suggests that in the linear regime of the CFI, for r S> 2200 ;^^^ the CFI spatially grows 
beyond the density ramp size L and spatially saturates with Lgat — Qt ^ L. 

Based on our analysis we further observe that the spatial-temporal nature of the instability also 
has an effect for finite beam-plasma interaction time Tint and beam size Abeam- la fact, depending 
on the relation between these parameters, i.e. either Abeam < QRnt or > a weaker or 
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stronger filamentation of the beam can be expected which will further affect the nonlinear 
growth of the instability. Moreover, the study of spatial-temporal evolution of the beam-plasma 
instabilities can also lead to a better understanding and characterization of the Weibel mediated 
collisionless shock formation process in laboratory and astrophysical plasmas. 
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